Hot carrier circuit reliability simulation

ABSTRACT

The present invention is directed to a number of improvements in methods for reliability simulations in aged circuits whose operation has been degraded through hot-carrier or other effects. A plurality of different circuit stress times can be simulated within a single run. Different aging criteria may be used for different circuit blocks, circuit block types, devices, device models and device types. The user may specify the degradation of selected circuit blocks, circuit block types, devices, device models and device types independently of the simulation. Device degradation can be characterized in tables. Continuous degradation levels can be quantized. Techniques are also described for representing the aged device in the netlist as the fresh device augmented with a plurality of independent current sources connected between its terminals to mimic the effects of aging in the device. The use of device model cards with age parameters is also described. To further improve the circuit reliability simulation, a gradual or multi-step aging is used instead of the standard one step aging process. Many of these features can be embedded within the circuit simulator. A user data interface is also presented to implement these techniques and further allow users to enter their device models not presented in the simulator. For example, a proprietary model of, say, the substrate current in an NMOS could used be with a SPICE simulator employing a different model to simulate the aging of the circuit.

CROSS REFERENCE TO RELATED APPLICATIONS

[0001] This application claims priority from provisional U.S. patent applications: serial No. 60/236,865, entitled “Hot-Carrier Circuit Reliability Simulation”, filed Sep. 29, 2000; 60/236,587, entitled “Hot-Carrier Device Degradation Modeling and Extraction Methodologies”, filed Sep. 29, 2000; and 60/236,586, entitled “Hot-Carrier Reliability Design Rule Checker”, filed Sep. 29, 2000.

[0002] The above referenced applications are incorporated herein by reference for all purposes. The prior applications, in some parts, may indicate earlier efforts at describing the invention or describing specific embodiments and examples. The present invention is, therefore, best understood as described herein.

FIELD OF THE INVENTION

[0003] This invention relates generally to a design aiding apparatus and methods for the design of integrated circuits, and, more specifically, to methods for simulating the time degradation of a circuit due to the hot carrier effect.

[0004] Copyright Notice

[0005] A portion of the disclosure of this patent document may contain materials that are subject to copyright protection. The copyright owner has no objection to the facsimile reproduction by anyone of the patent document or the patent disclosure, as it appears in the Patent and Trademark Office patent files or records, but otherwise reserves all copyright rights whatsoever.

BACKGROUND OF THE INVENTION

[0006] Complementary metal oxide semiconductor (CMOS) circuits may contain thousands or millions of transistors and other circuit elements. The design of such circuits is extremely complicated and usually employs simulators, such as SPICE, to predict circuit operation. In addition to predicting the operation of the circuit when it is new, it is frequently important to simulate the operation of the circuit as it ages from processes such as hot carrier effects.

[0007]FIG. 1 is a schematic of a field effect transistor, the example here being a P-channel MOSFET formed on a silicon substrate, although much of the following discussion applies NMOS transistors and other insulated gate devices. This device may also have a lightly doped drain (LDD) region that is not shown in FIG. 1. Some of the particulars for the NMOS case will be discussed later. The transistor 10 has a bulk N-type region 5, source region 11 and drain region 12 defined at and below the substrate surface, and a lightly doped N-type channel region 15 having a manufactured length L_(M) in the region between the source 11 and the drain 12. Over the channel region is an oxide or other insulating layer 14 that supports the gate 13. (Here the gate is shown doped as N+, although in a PMOS this may alternately be P+.) When the transistor is turned on, current flows through the channel region below the oxide/silicon interface.

[0008] With operation over time, the PMOS transistor 10 will suffer degradation due to device aging. An important mechanism in the aging of insulated gate devices is the “hot carrier” effect. When the transistor is turned on, a current I_(ds) of energized charge carriers flows through the channel 15 from the source 11 to the drain 12. Impact of these charge carriers with silicon atoms near the channel/drain juncture generates electron-hole pairs. Some of the resultant electrons have sufficient energy to pass through the insulating layer 14 by tunneling or other leakage mechanisms and are then collected at the gate 13, resulting in a gate current I_(G). The rest of the electrons result in a substrate current I_(sub) contribution, with generated holes adding to the drain current I_(ds). Some of the impact induced electrons become trapped inside of the insulating layer 14 or become trapped at the interface of the channel 15 and the insulating layer 14. This attracts holes to the channel side of the channel/insulator interface, which cause the effective boundary of the P+ drain region 12 to shift closer to the source region. This is shown in FIG. 2, where drain is extended by the region12′, shortening the effective channel length to L′<L_(M). Additionally, the oxide/silicon interface above channel region is damaged by the more energetic holes and electrons.

[0009] These effects degrade the transistor's operation in several ways. As these charges accumulate, the voltage needed at the gate to turn on the transistor, the device's threshold voltage V_(th), changes, resulting in positive shift, ΔV_(th)>0. This makes the PMOS threshold voltage less negative and reduces the V_(th) value at which the onset of leakage current is seen. When a sufficient number of trapped electrons are distributed over the manufactured length L_(M) of the channel 15, a relatively large leakage current will flow through the transistor even when the gate voltage V_(G) is set to an “off” level.

[0010] The damage to the interface results in a lower mobility, μ, as the carriers flow form source 11 to drain 12. The lower mobility results in more resistance and lowers I_(ds), where as the electrons trapped in the dielectric tend to raise the PMOS's I_(ds) curve. The strength of these effects depend differently on bias conditions, but traditionally the change due to lower mobility has been smaller; however, as device sizes have decreased, this interface damage becomes more significant and, depending upon bias conditions, often becomes the larger effect below a quarter micron or so.

[0011] These problems are aggravated as device sizes decrease into the submicron region. One reason is that an amount of incursion of the region12′ that produces a relatively small relative change for a channel length of, say, L_(M)≈1.5 μm becomes a much larger relative change at L_(M)≈0.25 μm. For example, if12′ extends an absolute distance of 0.05 μm, this produces a change of

(L′−L _(M))/L _(M)=−20%

[0012] for the shorter channel length, but of only −3% at the longer length. Another reason is that as the electric field in the channel is given by

E=V _(ds) /L′

[0013] the resultant field strength, and consequently the number of electron/hole pairs produced, increases greatly in a submicron device even at low power operation. For example, in a L_(M)≈0.25 μm MOSFET operating at 2.5 volts results in fields of ≈10⁷ V/m. AsL′ departs further and further from L_(M), the field strengths increase resulting in an even greater electron/hole pair production rate.

[0014] In the case of an NMOS transistor, the results of hot-carriers differ in several respects. Electrons will again be trapped in the dielectric and as the number of electrons in the oxide increases, a positive shift in voltage will again be needed at the gate to turn on the transistor results, ΔV_(th)>0. As V_(th) is positive now, the threshold voltage becomes larger in magnitude and the I_(ds) curve is lowered as the gate voltage must additionally overcome the negative charge embedded in the dielectric to cause the channel inversion.

[0015] The damage to the interface again lowers the mobility, μ, as the electrons flow form source 11 to drain 12, resulting in more resistance and lowering I_(ds). Thus in this case, both effects move the I_(ds) curve in the same direction. Also as with the PMOS, the strength of these effects depends differently on bias conditions, with the change due to lower mobility traditionally smaller but becoming increasingly important as device size shrink. Unlike the PMOS, though, the effective channel length does not decrease.

[0016] These hot electron effects build up over time causing the device to degrade as it ages. To determine the aging of a circuit as a whole, the degradation of the various devices within the circuit must be considered. As an example, consider the simple circuit shown in the schematic of FIG. 3, consisting of three inverters manufactured to be identical.

[0017]FIG. 3 shows three inverters, 20 a, 20 b, and 20 c, connected in parallel with a capacitive load C 25 a, 25 b, and 25 c connected to the output of each. Each of these inverters consists of a circuit like that of FIG. 4, with a PMOS transistor P_(a) 21 and an NMOS transistor N_(a) 22 connected as shown. To simulate the fresh circuit, a circuit simulation, such as of the SPICE or timing simulation type, can be performed using the same model card for each of PMOS transistors 21 and the same model card for each of the NMOS transistors.

[0018] To simulate the aged circuit, a model card representing the aged behavior of each of the elements needs to be used. To determine how the aged device will operate, the device is stressed to obtain the aging model information from the electrical test data. From this, the aged model card for the device can be extracted. However, the same aged model card can no longer be used for all of the similar devices. To see why, consider FIG. 3 again.

[0019] The amount of degradation in the PMOS 21 and the NMOS 22 of inverter 20 a depends on the input signal supplied at node X, while the degradation in inverters 20 b and 20 c will instead depend on the input at nodes Y and Z, respectively. The waveforms at these nodes will differ from each other as the input waveform will change in both shape and magnitude as it propagates through the circuit. Thus, if the operation of a circuit after, say, 5 years, of use is simulated, the relative amount of degradation in each of the inverters will differ and, consequently, a different aged model card will be needed for each of the transistors. For even the simple circuit of FIGS. 3 and 4, this results in a three-fold increase in the required number of model cards for the aged circuit simulation compared with the fresh circuit. As a real circuit will often contain thousands or even millions of elements, the simulation of the aging of such circuits is a extremely complicated process for which a large number of improvements are desirable.

SUMMARY OF THE INVENTION

[0020] The present invention is directed to a number of improvements in methods for reliability simulations in aged circuits whose operation has been degraded through hot-carrier or other effects. According to a first aspect, the circuit operation at a plurality of different circuit ages or stress times can be performed within a single simulation, greatly reducing the time required over what would be needed if each circuit stress time were to be simulated separately.

[0021] Another aspect of the present invention allows for different aging criteria to be used for different elements. For example, the circuit can be broken down into blocks, such as an analog sector and a digital sector which, due to their differing functions, would have their degradation judged by different aging criteria. Alternately, or additionally, different degradation standards or aging criteria can be selected for different devices, such as PMOS compared to NMOS, or device types, such as a MOSFET compared to diode or a bipolar junction transistor. Similarly, different criteria can be used for different device models, different circuit blocks, or different circuit block types.

[0022] Other aspects which produce an increase in simulation efficiency allow for the user to specify the degradation of selected circuit elements independently of the simulation and to store the values of intermediate parameters in a lookup table. The table may be pre-specified, produced in an earlier simulation, or be the result of a partial simulation ran for the purpose of determining table values.

[0023] Techniques are also described for representing the aged device in the netlist as the fresh device augmented with a plurality of independent current sources connected between its terminals to mimic the effects of aging in the device, where the strength and connections of the current sources are based on modeling different aging mechanisms. The continuous value device degradations can be quantized to build the aged device models, to build current sources, or both, so that rather that doing the simulation for the continuous aging, it can be done in a series of quantized steps. The use of device model cards with age parameters is also described.

[0024] To further improve the circuit reliability simulation, a gradual or multi-step aging is used instead of the standard one step aging process. This gradual aging techniques includes the effects due to the degradation of I_(ds), I_(sub), and other parameters in determining how the elements of the circuit age.

[0025] In a further aspect of the present invention, the device characterization means, the device estimation means, the device degradation quantization means, the circuit netlist generation for aging simulation, and other features can be embedded into the circuit simulator. This will improve simulation efficiency and capability still further.

[0026] A user data interface is also presented to implement these techniques and further allow users to enter their device models not presented in the simulator. For example, a proprietary model of, say, the substrate current in an NMOS could be used with a SPICE simulator employing a different model to simulate the aging of the circuit.

[0027] Additional aspects, features and advantages of the present invention are included in the following description of specific representative embodiments, which description should be taken in conjunction with the accompanying drawings.

BRIEF DESCRIPTION OF THE DRAWINGS

[0028]FIG. 1 is a schematic drawing of a MOSFET transistor.

[0029]FIG. 2 shows the effects of aging due to the “hot-carrier” effect in the device of FIG. 1.

[0030]FIG. 3 is a schematic of a simple circuit used to discuss the effects of aging on a circuit.

[0031]FIG. 4 is a detail of the circuit of FIG. 3.

[0032]FIG. 5 is a flow chart of hot-carrier circuit simulation.

[0033]FIG. 6 shows multiple current sources to represent different degradation mechanisms for a MOSFET.

[0034]FIG. 7 is a block diagram of an interface for use in hot-carrier or other aging simulation.

[0035]FIG. 8 shows the entry of MOSFET in a netlist converted to a parameterized aged device model subcircuit call.

DESCRIPTION OF REPRESENTATIVE EMBODIMENTS

[0036] The various aspects of the present invention will be described with respect to the flow chart of FIG. 5. As described in the following, this flow differs from the prior art methods in some of the steps included as well as in the details of some of the steps that are found in the prior art. Although the discussion is mainly in terms of aging due to hot carriers, the techniques generally apply to circuit degradation due to the degradation of its constituent devices from other mechanisms.

[0037] Step 101 provides the netlist describing the circuit that will have the simulation controlling commands and reliability parameters of the present invention. This will include the various connections of the circuit elements and model cards for these elements as well as the age models. The age models may be those that are described in provisional U.S. patent application serial No. 60/236,587, entitled “Hot-Carrier Device Degradation Modeling and Extraction Methodologies”, filed Sep. 29, 2000, which is hereby incorporated herein by this reference. In one improvement over the prior art, multiple hot-carrier circuit simulations can be accomplished in one run with one circuit simulation by setting proper hot-carrier simulation commands. This will improve simulation efficiency.

[0038] In the prior art flow, each run of a hot-carrier circuit simulation would require a complete run of all of the steps in the flow. For example, to simulate the operation of the circuit and determine its degradation after, say, 5 years and 10 years, a complete simulation would be performed for each of these values. Referring to FIG. 5, step 107 typically consumes by far the most time. By combining the simulation for multiple values, all of these stress times can be simulated with just one run of step 101-107. Although steps 109, 111, and 113 will each be performed for each parameter, by only requiring one simulation for the early steps—and particularly step 107—efficiency is greatly improved.

[0039] Another improvement of the present invention found in step 101 is that different device lifetime criteria can be applied for different circuit blocks, circuit block types, devices, device models, and device types in one or multiple hot-carrier circuit simulations. Examples of different lifetime criteria include using parameters such as I_(ds) at different bias conditions, device transconductance, threshold voltage, model parameters, and so on for determining the amount of degradation. These criteria are parameters measurable before and after a circuit element is stressed and that are computable to use as a gauge of device degradation. The different criteria can be important in different areas of the circuit. For example, a circuit may contain both an analog block and a digital block of elements. The important degradation characteristics of these blocks may differ, with, say, transconductance G_(m), V_(th), or some other voltage behavior being used in the analog elements while for a digital part it would be driving capability of I_(ds) or the circuits delay or speed. The different criteria can be important in determining individual device lifetime and setting a threshold value below which the device becomes unacceptable for circuit performance requirements, for example, when I_(ds) or G_(m) at some bias is degraded below a certain value.

[0040] Even when circuit blocks are of the same type, such as both being adders or inverters, different criteria can be applied if the degradation of the blocks can best be quantified by different parameter values. These criteria may be different device parameters, such as the gain G_(m) or threshold values V_(th), or the same parameter but with different values, say a V_(th) shift of 100 mV in one block and 200 mV in another.

[0041] Similarly, differing criteria can be used for different devices, with say one criterion for PMOS and one for NMOS as these may have differing degradation mechanisms; for different models of the same device, such as an NMOS with different model cards due to, say, different sizes; or for different device types, as the important criteria for a MOSFET differ from those of a diode, which both differ from those of a bipolar junction transistor, and so on. For example, while driving capability may be important for both NMOS and PMOS devices within a given circuit block, leakage current may only be important for a PMOS in another block; or while driving capability may be the important quantity in MOSFETs, in a bipolar junction transistor the important performance criterion may be, say, emitter-base leakage current or current gain, β=I_(c)/I_(d).

[0042] At step 103, the device degradation screening means allows a user to assign different amounts of device age, degradation and/or lifetime for different circuit blocks, circuit block types, devices, device models, and device types in absolute or relative scales. For those devices which have pre-specified degradations, the computation of their degradations can be skipped. Relevant device characterizations (shown in step 105) can also be skipped. This will improve simulation efficiency and reduce memory usage. This step is optional.

[0043] For example, in a circuit with an analog block and a digital block, it may be found that the analog block does not limit circuit life. The user could set this block as already having a specific absolute degradation value for better speed. Alternately, as a PMOS device usually suffers less degradation than an NMOS, the user could assign a smaller relative scale: If an NMOS ages amount x, the corresponding PMOS would age α·x , where α<1 is a proportionality parameter. In the circuit of FIGS. 3 and 4, setting a relative PMOS value would cut the transistor degradation calculations in half as the value used for PMOS 21 would be obtained from that of NMOS 22 instead of requiring a independent calculation.

[0044] The device degradation characterization means of step 105 is an optional step that characterizes the device degradation building tables which can be used during hot-carrier circuit simulation for device age, degradation and lifetime computations. The tables can be built during the hot-carrier circuit simulation or as part of a stand-alone process. In a stand-alone, a new table can be produced as part of a partial simulation ran for this purpose or be provided from another source. The tables can also be re-used for further simulations.

[0045] Different forms of table can be built. The ART (Age Rate Table) approach is one of the forms which builds a device age rate table under various bias conditions and which makes the device age calculation more efficient. Some other forms, such as an I_(sub) table, an I_(gate) table, and I_(sub)/I_(ds) table for MOSFET transistors eliminate the need of I_(sub) and/or I_(gate) calculation which reduces the computation cost of MOSFET device degradation.

[0046] For example, the change in device age over a time δt due to hot-electrons may be expressed in one model for an NMOS as ${{\delta \quad {{Age}(t)}} = {\frac{I_{ds}}{WH}\left( \frac{I_{sub}}{I_{ds}} \right)^{m}\delta \quad t}},$

[0047] where H and m are the model's device degradation parameters and W is the devices channel width, and for a PMOS as ${{\delta \quad {{Age}(t)}} = {\left\lbrack {{W_{g} \times \frac{1}{H_{g}}\left( \frac{I_{gate}}{W} \right)^{m\quad g}} + {\left( {1 - W_{g}} \right) \times \frac{I_{ds}}{WH}\left( \frac{I_{sub}}{I_{ds}} \right)^{m}}} \right\rbrack \delta \quad t}},$

[0048] where H_(g) and mg are degradation parameters to include the gate current of the PMOS and W_(g) is a weighting parameter. The age rate in the NMOS aging model above is then ${{{rate}\frac{\delta \quad {{Age}(t)}}{\delta \quad t}} = {{\frac{I_{ds}}{WH}\left( \frac{I_{sub}}{I_{ds}} \right)^{m}} \neq {f(t)}}},$

[0049] that, for a sufficiently small δt, is a time invariant quantity that can be saved as a table. These formulae show how the current ratios are related to the device's age and that the larger the stress time t, the greater the age. For example, as these currents are functions of the bias voltages, I_(ds), I_(sub)=f(V_(ds), V_(gs), V_(bs)), these intermediate parameter values can be saved in a three dimensional table of V_(ds), V_(gs) and V_(bs), or in a reduced dimensionality table if, say V_(bs) is less important and the user wants to save memory. This intermediate quantity table can then be used whenever the parameter comes up. This can be implemented as a command, say .table, to specify the values,

.table Vds_start, Vds_stop, Vds_step, Vgs_start, Vgs_stop, Vgs_step

[0050] or set up structure for the table. Tables can be in even or non-even steps.

[0051] The table can be provided before the simulation from an outside source, such as test data, since it is only for single device and does not need the netlist. Alternately, a simulation can be run just for producing such a table. The table size can be chosen by the user as this is a trade off between memory size and speed/accuracy. However, once the table is produced it can be saved for future use and not be rebuilt every time it is required. Additionally, the table can be built incrementally, with entries added for better resolution. The use of tables can also improve the simulation accuracy. If the table size is large enough, it will be more accurate than calculating on the fly.

[0052] The device degradation characterization means can be embedded into the circuit simulation of step 107. If a user has their own simulator for step 107 and source code or interface functions, step 105 can be combined into step 107. The advantage of keeping step 105 separate from step 107 is that the steps outside of 107 can be independent of simulator type.

[0053] Step 107 can be a standard prior art simulation and uses the fresh model. The possible circuit simulations could, for example, be a SPICE simulator (such as HSPICE or Spectre) or a timing simulator (such as Starsim or Timemill). As this produces the fresh waveform at all the nodes of the circuit, it determines the relative stress on the circuit elements. Memory usage can be minimized by virtue of the screening procedure as shown in step 103 and the characterization procedure as shown in step 105 as these reduce the number of waveforms that need to be stored. If step 107 stands alone, the process is more flexible as it can be used with the various prior simulators and the other steps can be connected through the interface as described below.

[0054] Using a separate simulator for step 107 is less efficient as it requires moving large amounts of data back and forth across the interface. If the simulation can be merged with other steps, it saves on both the time and storage required for this transferal. Embedding into the circuit simulator some, or all, of the functions from the device degradation characterization means (step 105), the device degradation estimation means (step 109), the device degradation quantization means (step 111), and the circuit netlist generation for aging simulation (step 113) will make the hot-carrier circuit simulation much more efficient and capable. This would let the whole flow to be performed as a single run, processing the data along with simulator and not requiring storage of all the other data. This can increase the total speed and capabilities of the flow by orders of magnitude. This embedded process can also combine the parameterized aged device model as described below. Note that this embedding can be done at the device level and not just at the block level.

[0055] Step 109 is the device degradation estimation to establish all of the ages for each transistor, or other aged device, specified in the netlist. The device degradation analysis will be more efficient because the relevant storage and computations for the devices pre-specified in step 103 can be skipped. Many devices suffer from multiple aging mechanisms that may have different relative importance depending upon bias conditions or other factors. For devices with multiple hot-carrier degradation taking effect, these devices will have multiple independent device age values each corresponding to one mechanism. These various age values will contribute to the total device degradation.

[0056] As already noted, the device degradation estimation means can be embedded into the circuit simulator. This may be implemented in the model evaluation module, but in this case the process depends upon the particular models employed. In one aspect of the present invention, the device degradation estimation means is embedded into the circuit simulator outside of the model part and is consequently independent of the particular models used. This allows newer device models to be incorporated without requiring the device age calculation having to embedded individually for each model.

[0057] Some examples of different aging mechanisms are the electron trapping and generation of interface state at the oxide/silicon interface. In addition to trapping electrons, the oxide may also capture holes in traps. Due to their differing charge, trapped holes will react to the resultant electric field at a bias level differently from the electrons: If the drain to gate voltage is positive, V_(dg)=V_(d)−V_(g)>0, the holes will migrate towards the gate and the electrons towards the drain. Thus, only holes can jump the interface at this bias, while the situation is reversed if V_(dg)<0.

[0058] At step 111, device degradation can be represented in the new circuit netlist by using the aged device models, by adding current sources, or both. In the first case, the elements and connections in the SPICE (or other simulation) netlist are not actually changed, but the aged model cards are used. As the voltage waveforms supplied to different elements are not the same, elements that have the same model card in the fresh simulation will generally have different aged model cards as the parameters will now differ. In the second case, connections of the netlist are changed by introducing current sources between device nodes, with independent current sources accounting for different degradation mechanisms. Examples of models which can be used for these current sources for MOSFETs are the DeltaMOS, DeltaLogId, and DeltaId models described in provisional U.S. patent application serial No. 60/236,587, entitled “Hot-Carrier Device Degradation Modeling and Extraction Methodologies”, which was included by reference above.

[0059] Adding multiple current sources to represent different degradation mechanisms is shown in FIG. 6 for a MOSFET with multiple sources used simultaneously. In this example, for MOSFET 60 which has multiple hot-carrier degradation mechanisms, each mechanism will have a corresponding device age (as shown in step 109), and consequently, a current source I_(i) 64-i, 1≦i≦n, will be added in parallel between the source S 62 and drain D 63. In this case, the aged device is represented in the new netlist by the fresh device between the nodes S′ 62′ and D′ 63′ and the aging effects separated out in the currents running from between the original source connection S 62 and the node S′ 62′ to between the original drain connection D 63 and the node D′ 63′. Each of the multiple current sources,

I _(i) =f _(i)(V _(ds) ,V _(gs) ,V _(bs) ;P),

[0060] is added to cover the different mechanisms, where the functional form depends on the mechanism and where P represents the parameter set for the model, such as that described in provisional U.S. patent application serial No. 60/236,587, entitled “Hot-Carrier Device Degradation Modeling and Extraction Methodologies”. The functional forms of the various f_(i) will be from physically based models with parameters set by electrical test data.

[0061] In some case, a single mechanism may be represented by more than one current source. For example, although the channel in a fresh transistor may be symmetric between the source and drain, due to damage at the silicon-oxide interface this symmetry is lost as the circuit ages. Consequently, the current sources added for the aged device may include separate sources for forward and reverse biased conditions.

[0062] The continuous value device degradations can be quantized to build the aged device models, to build current sources, or both, so that rather than doing the simulation for the continuous aging, it can be done in a series of quantized steps. For example, say a circuit has a million transistors that have the same model card when new. As each of these will generally have a different age ranging from zero to a maximum value, the aged device simulation would have a million aged cards in the continuous case. In the quantized case, the age range from zero to the maximum value is broken up into, say, ten sub-ranges with a single age value used for each. Then only the discrete subset of ten model cards is used. The full interval need not be broken up into sub-ranges of uniform size if better resolution is more important in a particular range. The criteria for model quantization can be device age or any device parameters (such as I_(ds) at different bias conditions, device trans-conductance, threshold voltage, etc.) degradations. The device degradation quantization means can be embedded into the circuit simulator. This step is optional.

[0063] In step 113 the aging circuit simulation is run with the different elements at their respective different aged models. Each of the simulations of step 113 will generally take significantly less time than the simulation already performed for the fresh device in step 107. This is because in running step 107, the simulation will record the waveforms at the required nodes. These nodes can be specified for particular devices in the screening stage. In step 113, only waveforms at the nodes that are needed for comparison with the fresh waveforms need to be stored. Thus, when doing a simulation for multiple ages in a single run, the ability to share steps 101-107 results in large increase in efficiency even though steps 109-113 must be performed individually for each age.

[0064] The above flow is the so-called “one step aging” circuit reliability simulation, where device ages and degradation are estimated using the fresh waveforms generated in step 107. However, as the circuit ages, the waveforms at the devices comprising the circuit will change. Therefore, the accuracy of the circuit reliability simulation should be improved if the circuit waveforms are updated during the degradation circuit simulation. This process will be called “gradual aging”.

[0065] “Gradual aging” can be realized by adding a loop for steps 107, 109 and 111, as shown by the dotted line in FIG. 5, to update the models of the circuit simulator. Assuming the total circuit stress time is T_(n), e.g., 10 years, T_(n) is divided into N time intervals, [T₀, T₁], [T₁, T₂], . . . ,[T_(n−1), T_(n)], where T₀ is zero. The transient simulation time is T, which is usually much smaller than any above time intervals, for example 10 μs. In “one step aging”, the circuit waveform from [0, T] is projected to the end of T_(n), ${{{Age}\left( T_{n} \right)} = {\frac{T_{n}}{T}{{Age}(T)}}},$

[0066] or, for the 10 year, 10 microsecond example, ${{Age}\left( {10\quad {yr}} \right)} = {\frac{10\quad {yr}}{1\quad \mu \quad s}{{{Age}\left( {1\quad \mu \quad s} \right)}.}}$

[0067] Instead, for gradual aging the projection will now be made for each time interval.

[0068] At the beginning of each time interval, device models will be updated according to device age values. Circuit simulation 107 will be launched using updated device models in each loop, with only the first time interval using the fresh device model. Device age and degradation will then be estimated for each loop in step 109. Device degradation can optionally be quantized and a new circuit netlist for the next time interval will be prepared in step 111. Thus, to make the “gradual aging” more accurate for MOSFET devices, in addition to making I_(ds)=I_(ds)(t), degraded models for I_(sub), I_(gate), and intermediate values can be used for device degradation calculation. For example, in the NMOS hot-carrier aging model example of step 105, age is now a function of I_(sub)(t) and I_(ds)(t): ${{Age}(t)} = {\sum\limits_{i}{\frac{I_{ds}\left( T_{i} \right)}{WH}\left( \frac{I_{sub}\left( T_{i} \right)}{I_{ds}\left( T_{i} \right)} \right)^{m}\Delta \quad {t_{i}.}}}$

[0069] This loop will continue until the final time interval is completed. While the flow calculated the device age parameters (H, m) using I_(ds), I_(sub), and so on at step 107 for the fresh simulation, the loop back from step 113 uses a differing model at step 107 as I_(ds), I_(sub), and the other parameters will change with time.

[0070] Open Data Interface

[0071]FIG. 7 is a block diagram of an open data interface (ODI) for use in hot-carrier or other aging simulations. The interface 203 is connected between the degradation simulator 201 and the user inputs 205. In block 205 some of information the user may supply is shown, although not all of these will necessarily be provided in all cases, and in some instance will be integrated into the simulator as described above with respect to step 107 of FIG. 5. This interface 203 places the hot-carrier circuit simulator 201 into the designer's environment, allowing the user to control the calculation of age and other functions by supplying their own, say, age equations and other models into the simulator through 205 to calculate age and degradation equations. This allows the user to customize the simulation to their needs.

[0072] Referring to block 205 in more detail, user defined data and models (such as device age 209, device degradation 211 and auxiliary functions 207) that will be used for device age and degradation calculation and user defined circuit simulator can communicate with the hot-carrier circuit simulator through ODI (203). The Circuit Simulator 213 can be of the SPICE or timing simulator type, or the user's adaptation or equivalent. Device instance and model parameters are passed from the hot-carrier circuit simulator 201 through the ODI 203 to user defined data 205. For example, if the channel width W of a MOSFET is needed for the age calculation or the bias conditions are needed for the auxiliary functions, the hot-carrier circuit simulator 201 can provide these values. Similarly, user defined data and the simulation results from the user defined circuit simulator are passed back to hot-carrier circuit simulator 201 through ODI 203.

[0073] The user can also supply auxiliary functions to replace or supplement the equations found in the hot-carrier circuit simulator 201. For example, I_(gate) and I_(sub) are sample auxiliary functions. The simulator will provide default equations, but the user may provide auxiliary function that are lacking or for which the user has a proprietary model that they prefer over the default model for their application. For the example of I_(sub), many different model for the substrate current exist and the user may prefer a particular one in their application, such as described in co-pending U.S. patent application No. 09/661,328, filed on Sep. 14, 2000, entitled “MOSFET Modeling for IC Design Accurate for High Frequencies”, which is hereby incorporated herein by this reference. The user can select which function they supply, just specifying a module for, say, I_(sub), and using default modules for others.

[0074] Input parameters, intermediate and final results can be shared among all the user defined data, such as passing I_(sub) to I_(g), as I_(g)∝I_(sub). This can also include bias conditions, so that a function, say, age can now be a function of parameters beyond I_(ds), age→age(I_(ds),V_(ds),V_(gs),V_(bs)), and so on. Multiple user defined data of the same functionality (such as the device age) can be integrated with the other user defined data. The user may have multiple equations of how, say, age model is calculated and may call any of these. Users have the flexibility to choose any one of them. User defined data can also be combined with the internal data in the hot-carrier circuit simulator (201) to fulfill the device age/degradation calculation. Again, these can be data that come with device or add in some proprietary ones the user has added in. The user defined data can be incomplete. For any of the user definable data not defined, internal data will be used. User defined functions can access intermediate results of circuit simulation, such as the saturation voltage, V_(dsat), or threshold voltage, V_(th), and not just I_(ds).

[0075] Parameterized Aged Device Model

[0076] Another approach that can be used to study circuit degradation is the use of a parameterized device model. This method can be used to replace the above approach through use of a parameterized model card for each of the specified circuit elements, or can be used in a complementary fashion, with some of devices using the parameterized model as described below and the rest of the devices treated as above. Each model parameter can be formulated as a function of several intermediate quantities that measure device degradation, such as device age or other device parameter (I_(ds), V_(th), G_(m), and so on). We will use device age as an example hereinafter.

[0077] Because aged device parameter is a function of device age, device model card can be parameterized with the device age as a parameter. This implies that one model card can be used by devices of different device age values. Age value can be quantized to limit the number of aged device models if a circuit simulator needs to predetermine the models at the start of simulation. One sample model card could be expressed as

.model nmos nmos vth0=‘0.7+5*age’,

[0078] where vth0 is the model parameter related to the threshold value and the other model parameters are suppressed for simplicity. The functional form of the parameters, such as the shown linear relation for the threshold voltage, can either be experimentally determined, or have a functional form based on physical model with the coefficients experimentally determined. The building of age models is described more fully in provisional U.S. patent application serial No. 60/236,587, entitled “Hot-Carrier Device Degradation Modeling and Extraction Methodologies”, which was included by reference above.

[0079] The circuit simulator can then accept device age as an instance parameter. One example is listed as follows:

m1 1 2 3 4 nmos W=10u L=0.25u age=0.01

[0080] This age value will be used for calculating aged device model parameters. To obtain the age values for the devices in the netlist, the steps 101-109 in the flow of FIG. 5 can still be performed, but the preparation of the netlist in step 111 is different. The age value can also be specified.

[0081] In simulators that do not support age as an instance parameter, a macro or subcircuit can also be used to prepare aged device model parameters from a device age value. An auxiliary program can convert each degradable device to a subcircuit call. The advantage of this approach is that currently available circuit simulators widely support this netlist style even if they do not accept device age directly as an instance parameter. The examples of FIG. 8 shows the entry of MOSFET m1 in a netlist converted to a subcircuit call Xm1 with the device age value of age1. The left and right show two differing results, on the left converting to a subcircuit to introduce the age and on the right to compute a parameter (here vth0 as above) from this value. The value of age1 will be used to prepare the aged device model for m1. The use of subcircuits is described more fully in co-pending U.S. patent application No. 09/661,328, filed on Sep. 14, 2000, entitled “MOSFET Modeling for IC Design Accurate for High Frequencies”, that was incorporated by reference above.

[0082] According to another aspect, device age calculation and aged device model parameters can be accomplished inside the model evaluation module for each device. A circuit simulator contains both an engine to solve the circuit and the model evaluation module. Rather than compute an expression such as, say, vth0=‘0.7+5*age’ in the engine, the model evaluation module can alternately calculate vth0=f(age). A model flag can be used to control whether the computation of the device age and the aged device model parameter should be carried out. For example, this can be implemented

.model nmos nmos aging_flag=1,

[0083] where the other parameters are again suppressed. If the flag is set to 1, the aged device model parameter will be calculated based on the parameterized equation. This would occur in the circuit simulation 107 and present a more integrated approach to compute I_(ds) and other output values. If the flag is set to 0, the calculation would not occur, as would be the case in simulating the fresh model or when the user prefers to use the full flow of FIG. 5.

[0084] Many aspects of the methods of the present invention will most commonly be implemented in software as a computer program product, although many of these can be implemented in hardware or by a combination of software and hardware. As will be understood in the art, the invention or components thereof may be embodied in a fixed media program component containing logic instructions and/or data that when loaded into an appropriately configured computing device cause that device to perform according to the invention. As will be understood in the art, a fixed media program may be delivered to a user on a fixed media for loading in a users computer or a fixed media program can reside on a remote server that a user accesses through a communication medium in order to download a program component. Examples of such fixed media include a disk-type optical or magnetic media, magnetic tape, solid state memory, etc.. The invention may be embodied in whole or in part as software recorded on this fixed media.

[0085] The invention also may be embodied in whole or in part within the circuitry of an application specific integrated circuit (ASIC) or a programmable logic device (PLD). In such a case, the invention may be embodied in a computer understandable descriptor language which may be used to create an ASIC or PLD that operates as herein described.

[0086] Although the various aspects of the present invention have been described with respect to specific exemplary embodiments, it will be understood that the invention is entitled to protection within the full scope of the appended claims. 

It is claimed:
 1. A method of simulating the degradation of a circuit, comprising: providing a netlist specifying the components of the circuit; supplying a plurality of circuit stress time values; supplying aging model information on selected ones of the components; simulating the behavior of the fresh circuit to determine for each of the selected components a component degradation parameter relative to circuit stress time; and determining the degraded operation of the circuit by simulating the operation of the circuit with the specified components using their respective aging model information and respective relative component degradation parameter at the supplied circuit stress time values.
 2. The method of claim 1, wherein said degradation of the circuit is due to hot-carrier effects.
 3. The method of claim 1, wherein the simulating is performed using a SPICE type circuit simulator.
 4. The method of claim 1, wherein the simulating is performed using a timing simulation type circuit simulator.
 5. The method of claim 1, wherein the aging model information on the selected ones of the components is derived from electrical test data.
 6. The method of claim 1, wherein said simulating the behavior of the fresh circuit determines the waveforms at the nodes to which the selected ones of the components are connected relative to an input waveform.
 7. The method of claim 1, wherein determining the degraded operation of the circuit comprises determining the circuit's speed at the supplied circuit age parameters.
 8. A method of simulating the degradation of a circuit, comprising: providing a netlist specifying the components of the circuit, wherein the circuit includes a plurality of distinct sets of components; supplying an independent performance criterion for each set of said plurality of distinct sets of components; supplying a circuit stress time value; supplying aging model information on selected components from each of said sets of components; simulating the behavior of the fresh circuit to determine for each of the selected components a component degradation parameter relative to circuit stress time, wherein each of the selected components' relative degradation parameter is determined using the respective performance criteria of the set to which the selected component belongs; and determining the degraded operation of the circuit by simulating the operation of the circuit with each of the specified components using the respective aging model information and respective relative component degradation parameter at the supplied circuit stress time value.
 9. The method of claim 8, wherein the distinct sets of components each form different functional blocks.
 10. The method of claim 8, wherein a first of said sets of components is an analog block and a second of said sets of components is a digital block.
 11. The method of claim 10, wherein the performance criterion of the first set is transconductance and the performance criterion of the second set is drain to source current.
 12. The method of claim 8, wherein the distinct sets of components consist of different device types.
 13. The method of claim 12, wherein a first of said different device types is an NMOS and a second of said different device types is a PMOS.
 14. The method of claim 13, wherein the PMOS performance criterion is leakage current.
 15. The method of claim 13, wherein the NMOS performance criterion is driving capability.
 16. The method of claim 12, wherein a first of said different device types is a MOSFET and a second of said different device types is a bipolar junction transistor.
 17. The method of claim 16, wherein the bipolar junction transistor performance criterion is leakage current.
 18. The method of claim 8, wherein the distinct sets of components employ different models for simulating the same device type.
 19. The method of claim 8, wherein the distinct sets of components consist of the same device type.
 20. The method of claim 8, wherein the distinct sets of components form functional blocks performing the same function.
 21. A method of simulating the degradation of a circuit, comprising: providing a netlist specifying the components of the circuit; supplying a circuit stress time value; supplying aging model information on a first set of selected components of the circuit; simulating the behavior of the fresh circuit to determine for each of the first selected set of components a component degradation parameter relative to circuit stress time; specifying the degradation level of a second set of selected components of the circuit, wherein the elements of the first set and the second set of components are distinct; and determining the degraded operation of the circuit by simulating the operation of the circuit with each of the first set of specified components using the respective aging model information and respective relative component degradation parameter at the supplied circuit stress time value and with each of the second set of specified components using the respective specified degradation level.
 22. The method of claim 21, wherein the first and second sets of components each form different functional blocks.
 23. The method of claim 21, wherein the second set of components form a digital block.
 24. The method of claim 23, wherein the specified degradation level is expressed in terms of drain to source current degradation.
 25. The method of claim 21, wherein the second set of components is an analog block.
 26. The method of claim 25, wherein the specified degradation level is expressed in terms transconductance degradation
 27. The method of claim 21, wherein the first and second sets of components each consist of different device types.
 28. The method of claim 27, wherein the second set of components consists of PMOS transistors.
 29. The method of claim 28, wherein the specified degradation level is expressed in terms of leakage current degradation.
 30. The method of claim 27, wherein the second set of components consists of NMOS transistors.
 31. The method of claim 30, wherein the specified degradation level is expressed in terms of driving capability degradation.
 32. The method of claim 27, wherein the second set of components consists of bipolar junction transistors.
 33. The method of claim 32, wherein the specified degradation level is expressed in terms of leakage current degradation.
 34. The method of claim 21, wherein the first and second sets of components employ different models for simulating the same device type.
 35. The method of claim 21, wherein the first and second sets of components each consist of the same device type.
 36. The method of claim 21, wherein the first and second sets of components form functional blocks performing the same function.
 37. The method of claim 21, wherein the degradation level of the second set of selected components is specified as a relative component degradation parameter with respect to the component degradation parameter of the first set of components.
 38. The method of claim 21, wherein the degradation level of the second set of selected components is expressed in terms of age.
 39. The method of claim 21, wherein the degradation level of the second set of selected components is expressed in terms of lifetime.
 40. A method of simulating the degradation of a circuit, comprising: providing a netlist specifying the components of the circuit; supplying a plurality of circuit stress time values; providing a device degradation table; determining for selected ones of the components a component degradation parameter relative to circuit stress time from the device degradation table; and determining the degraded operation of the circuit by simulating the operation of the circuit with the specified components using their respective relative component degradation parameter at the supplied circuit stress time values.
 41. The method of claim 40, wherein the device degradation table is expressed in terms of one or more intermediate quantities.
 42. The method of claim 41, wherein the device degradation table is time independent.
 43. The method of claim 42, wherein the device degradation table is an age rate table.
 44. The method of claim 41, wherein the device degradation table is a function of bias voltages.
 45. The method of claim 40, wherein said providing a device degradation table comprises: supplying aging model information on selected ones of the components; simulating the behavior of the fresh circuit to determine for each of the selected components a component degradation parameter relative to circuit age; and storing an intermediate quantity derived from the relative component degradation parameters as said device degradation table.
 46. A method of simulating the degradation of a circuit, comprising: providing a netlist specifying the components of the circuit; supplying a circuit stress time value; providing a device degradation table expressed in terms of device current; determining for selected ones of the components a component degradation parameter relative to circuit stress time from the device degradation table; and determining the degraded operation of the circuit by simulating the operation of the circuit with the specified components using their respective relative component degradation parameter at the supplied circuit stress time value.
 47. The method of claim 46, wherein the device degradation table is expressed in terms of gate current.
 48. The method of claim 46 wherein the device degradation table is expressed in terms of substrate current.
 49. The method of claim 46, wherein the device degradation table is expressed in terms of the ratio of substrate current to drain to source current.
 50. A method of simulating the degradation of a circuit, comprising: providing a netlist specifying the components of the circuit; supplying a circuit stress time value; providing a device degradation table; determining for selected ones of the components a component degradation parameter relative to circuit stress time from the device degradation table; and determining the degraded operation of the circuit by simulating the operation of the circuit with the specified components using their respective relative component degradation parameter at the supplied circuit stress time value, wherein said simulating the operation of the circuit is performed with a circuit simulator and wherein said providing a device degradation table is embedded in the circuit simulator.
 51. A method of simulating the degradation of a circuit, comprising: providing a netlist specifying the components of the circuit; supplying a circuit stress time value; supplying aging model information on selected ones of the components; simulating the behavior of the fresh circuit to determine for each of the selected components a component degradation parameter relative to circuit stress time; revising the netlist, wherein each of said selected components is replaced by a non-aged version of the selected component and a plurality of independent current sources corresponding to different mechanisms connected between the terminals of the non-aged version, wherein the magnitude of the current relative to a circuit stress time in each of the current sources of a component is determined from the aging model information of the component and a distinct mechanism degradation parameter derived from the component degradation parameter; and determining the degraded operation of the circuit by simulating the operation of the circuit with the revised netlist at the supplied circuit stress time value.
 52. The method of claim 51, wherein the selected components are MOSFETs.
 53. The method of claim 52, wherein said degradation of the circuit is due to hot carrier effects.
 54. The method of claim 52, wherein for each of said selected components more than one of said plurality of independent current sources are connected between the source and drain terminals of the non-aged version.
 55. The method of claim 52, wherein said method further includes: determining the magnitude of the respective current in each of the independent current sources, said determining comprising: supplying a physical model of the current magnitude; and establishing the values of the coefficients in the physical model from electrical test data.
 56. The method of claim 52, wherein the degradation level of the selected components is expressed in terms of lifetime.
 57. The method of claim 52, wherein the degradation level of the selected components is expressed in terms of age.
 58. The method of claim 51, wherein said simulating the operation of the circuit is performed with a circuit simulator and wherein said revising the netlist is embedded in the circuit simulator.
 59. A method of simulating the degradation of a circuit, comprising: providing a netlist specifying the components of the circuit; supplying a circuit stress time value; supplying aging model information on selected ones of the components; simulating the behavior of the circuit to determine for each of the selected components a component degradation parameter relative to circuit stress time, wherein said simulating the behavior is performed using a circuit simulator and includes incorporating the aging of the selected components by updating the models of said circuit simulator; revising the netlist, wherein each of said selected components is replaced by a non-aged version of the selected component and a plurality of independent current sources corresponding to different mechanisms connected between the terminals of the non-aged version, wherein the magnitude of the current relative to a circuit stress time in each of the current sources of a component is determined from the aging model information of the component and a distinct mechanism degradation parameter derived from the component degradation parameter; and determining the degraded operation of the circuit by simulating the operation of the circuit with the revised netlist at the supplied circuit stress time value.
 60. A method of simulating the degradation of a circuit, comprising: providing a netlist specifying the components of the circuit; supplying a circuit stress time value; supplying aging model information on selected ones of the components; simulating the behavior of the circuit to determine for each of the selected components a component degradation parameter relative to circuit stress time, wherein said simulating the behavior is performed using a circuit simulator and includes incorporating the aging of the selected components by updating the models of said circuit simulator; and determining the degraded operation of the circuit at the supplied stress time value by simulating the operation of the circuit with the each of the specified components using the respective aging model information and respective relative component degradation parameter at the supplied stress time value.
 61. The method of claim 60, wherein the selected components are MOSFETs and said incorporating the aging of the selected components comprises including the time dependence of the drain to source current.
 62. The method of claim 61, wherein said incorporating the aging of the selected components comprises including the time dependence of the substrate current.
 63. The method of claim 61, wherein said incorporating the aging of the selected components comprises including the time dependence of the gate current.
 64. The method of claim 60, wherein said simulating the behavior of the circuit to determine for each of the selected components a component degradation parameter relative to circuit age comprises: simulating the behavior of the fresh circuit to determine for each of the selected components an intermediate component degradation parameter relative to circuit stress time; determining the degraded operation of the circuit at an intermediate circuit stress time value by simulating the operation of the circuit with the each of the specified components using the respective aging model information and respective relative intermediate component degradation parameter at the intermediate circuit stress time value; and simulating the behavior of the degraded circuit at the intermediate circuit stress time value to determine for each of the selected components a component degradation parameter relative to circuit age, wherein the intermediate circuit stress time value is less than the circuit stress time value.
 65. The method of claim 60, wherein said simulating the operation of the circuit is performed with a circuit simulator and wherein said incorporating the aging of the selected components by updating the models of said circuit simulator is embedded in the circuit simulator.
 66. A method of simulating the degradation of a circuit, comprising: providing a netlist specifying the components of the circuit; providing for selected ones of the components model cards containing a device degradation parameter; supplying a circuit stress time value; simulating the behavior of the fresh circuit to determine for each of the selected components a device degradation parameter value relative to stress time; and determining the degraded operation of the circuit by simulating the operation of the circuit with each of the selected components' model cards using the respective device degradation parameter corresponding to the supplied circuit stress time value.
 67. The method of claim 66, wherein said simulating the operation of the circuit is performed by a circuit simulator comprising an engine and a model evaluation module, and wherein said component degradation parameter is determined in the model evaluation module.
 68. The method of claim 66, wherein the device degradation parameter is an intermediate quantity.
 69. The method of claim 66, wherein the device degradation parameter is device lifetime.
 70. The method of claim 66, wherein the device degradation parameter is device age.
 71. The method of claim 66, wherein the device degradation parameter is device degradation.
 72. The method of claim 66, wherein the selected components model cards are implemented as a subcircuit.
 73. The method of claim 66, wherein the model cards containing a device degradation parameter contain one or more device model parameters as functions of said device degradation parameter.
 74. The method of claim 73, wherein the functions of said device degradation parameter are extracted from experimental data.
 75. A method of simulating the degradation of a circuit, comprising: providing a netlist specifying the components of the circuit; supplying a circuit stress time value; supplying aging model information on selected ones of the components; simulating the behavior of the fresh circuit to determine for each of the selected components a degradation level relative to circuit stress time; quantizing each of said relative degradation levels to one of a plurality of discrete values; and determining the degraded operation of the circuit by simulating the operation of the circuit with the specified components using their respective aging model information and respective quantized relative degradation level at the supplied circuit stress time value.
 76. The method of claim 75, wherein said determining comprises: revising the netlist, wherein each of said selected components is replaced by a non-aged version of the selected component and a plurality of independent current sources corresponding to different mechanisms with distinct quantized relative degradation level connected between the terminals of the non-aged version, the magnitude of the respective quantized current in each of the current sources determined from the aging model information of component; and determining the degraded operation of the circuit by simulating the operation of the circuit with the revised netlist, the independent current magnitudes derived form the respective aging model information and respective relative degradation level at the supplied circuit stress time value.
 77. The method of claim 75, wherein said simulating the operation of the circuit is performed with a circuit simulator and wherein said quantizing is embedded in the circuit simulator.
 78. A method of simulating the degradation of a circuit, comprising: providing a netlist specifying the circuit at the device level; supplying a circuit stress time value; supplying aging model information on selected ones of the devices of which the circuit is comprised; simulating the behavior of the fresh circuit with a circuit simulator; determining for each of the selected devices a degradation parameter relative to circuit stress time, wherein said determining is embedded in said circuit simulator and is model independent; and simulating the degraded operation of the circuit with the specified components using their respective aging model information and respective relative device degradation parameter at the supplied circuit stress time value.
 79. The method of claim 78, wherein said selected devices are MOSFETs.
 80. The method of claim 78, further comprising: quantizing each of said device degradation parameters to one of a plurality of discrete values, wherein said simulating the degraded operation is performed using the quantized device degradation parameters.
 81. The method of claim 80, wherein said quantizing is embedded in the circuit simulator.
 82. A system for a user to simulate the degraded operation of a circuit, comprising: a circuit degradation simulator; a user defined circuit simulator, wherein said user defined circuit simulator allows the user to input a user selected circuit simulator and user data comprising: a netlist representation of the circuit; and a circuit stress time value; and an interface for connecting the user defined circuit simulator to the circuit degradation simulator, whereby the circuit degradation simulator determines the operation of the circuit at the circuit stress time value utilizing said netlist representation in the user selected circuit simulator.
 83. The system of claim 82, wherein said user selected circuit simulator is a SPICE type simulator.
 84. The system of claim 82, wherein said user selected circuit simulator is a timing simulator.
 85. The system of claim 82, wherein said user data further comprises: a user selected device age model, wherein the circuit degradation simulator utilizes said user selected auxiliary function in determining the operation of the circuit at the circuit stress time value.
 86. The system of claim 82, wherein said user data further comprises: a user device degradation model, wherein the circuit degradation simulator utilizes said user device degradation in determining the operation of the circuit at the circuit stress time value.
 87. The system of claim 86, wherein said user data further comprises: a user selected auxiliary function, wherein the circuit degradation simulator utilizes said user selected auxiliary function in determining the operation of the circuit at the circuit stress time value.
 88. The system of claim 87, wherein intermediate results in determining the operation of the circuit at the device age parameter are exchanged between said user device degradation model and said user selected auxiliary function.
 89. The system of claim 88, wherein said user selected auxiliary function is the substrate current.
 90. The system of claim 88, wherein said user selected auxiliary function is the gate current.
 91. A computer readable storage device embodying a program of instructions executable by a computer to perform the method of any one of claims 1, 8, 21, 46, 50, 51, 59, 60, 66, 75, and
 78. 92. A method for transmitting a program of instructions executable by a computer to perform a process of simulating the degradation of a circuit, said method comprising: causing the transmission to a client device a program of instructions, thereby enabling the client device to perform, by means of such program, the process of the method of any one of claims 1, 8, 21, 46, 50, 51, 59, 60, 66, 75, and
 78. 